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ABSTRACT 

First- and second-order explicit difference schemes are derived for a 
three-dimensional hyperbolic system of conservation laws, without recourse to 
dimensional factorization. All schemes are upwind biased and optimally 
stable. 
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1. INTRODUCTION 


For the solution of initial-value problems governed by hyperbolic conservation 
equations some fine numerical techniques are available. Most methods are based on a 
one-dimensional scheme 



The convenience of such a factorization is twofold. Firstly, in developing 
multi-dimensional methods one may concentrate on one-dimensional operators; 
secondly, multi-dimensional codes reduce to a sequence of one-dimensional sweeps. 
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A display of the power of the above approach is found in the review paper by 
Woodward and Colella [ 2 ] on the numerical simulation of two-dimensional, strongly 
compressible flow. 


Considering the success of dimensional factorization we may ask ourselves if there 
is any point in designing genuinely multi-dimensional methods, i.e. methods that 
can not be implemented as a sequence of one-dimensional operators. The answer still 

is yes , but may very well tend to "no" if computers will continue to grow, in 

speed and capacity, at the current pace. The same, however, may have been said, 
more than 20 years ago, on the matter of developing second-order methods. In both 
cases the increase in complexity of the methods is meant to pay off via an increase 
in efficiency. 

Genuinely multi-dimensional schemes are most efficient in modeling essentially 

dimensional phenomena in a relatively small region, such as push-pull flow or 

flow around a sharp corner. If such regions have a strong Influence on the overall 

solution, then the use of a multi-dimensional scheme is recommendable. 

In the further sections of this paper I shall indicate how to construct second- 
order difference schemes in two or three space dimensions. First-order schemes were 
published earlier [ 3 ]; their construction is briefly reviewed. 


2. A FIRST-ORDER SCHEME 

In order to derive a first-order scheme for Eq.(3), consider the scalar linear 
equation 

u^ + au^ + bu^ + cu^ =0,a>0,b>0,c>0, ( 5 ) 

in combination with a piecewise uniform initial-value distribution 


n n 

u (x,y,z) = 


X, , < X < 

1 -^ 


‘i+v y " ^J+%- 


< z < z. 


( 6 ) 


with = x.±^,Ax , y^^j^ = y^±^Ay , z^^ = z^^i^Az. 

The exact solution of (5) at time t in terms of the initial values at t*^ is 
n+1, n, 

u (x,y,z) = u (x-aAt,y-bAt,z-cAt). ( 7 ) 

Inserting the initial values (6) into (7) yields, for sufficiently small At, the 
following difference scheme: 


"" ■ h - - T - i Yy>l 

- <,4 {1 - - i - W 4 (1 - i „ 4 J) 
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-aA{l-^aA(l-iaA)-^aA(l--YaA )}]u”, 

•7 7 ^ X X i \7 \7 V V A V 'X’ J J ' 


3 \j «-* 

y y 


y y 


3 v/ / 

X X 


a = aAt/Ax < 1 , a = bAt/Ay < 1 , a = cAt/Az < 1 ; 
X y z 


( 8 ) 


here A^, A^, A^ denote backward (upwind) differencing in the x-, y-, z-direction. 


The terms in the above formula have been especially arranged in order to facilitate 
their interpretation. The terms in parentheses are one-dimensional updates over a 
time-step 1/3 At, with the first-order upwind scheme. The braced terms are two- 
dimensional updates over a time-step \ At, and represent the values of u at cell 
interfaces, averaged both in space and in time. The full operator, between 
brackets, shows the three-dimensional upwind nature of the scheme. 

Godunov [ 4 ] has indicated how to transform upwind differencing for a scalar linear 
convection equation into a method for a nonlinear hyperbolic system of the form 
(2). For scheme (8) his recipe implies replacement of a A u by 


f (V>'1 


or 




(9.1) 

(9.2) 


where V denotes forward differencing and the superscripts + and - indicate the 
splitting of a flux-difference or a flux in parts associated with forward and 
backward signals, respectively (see [5] and [6]). The expressions <^yA u and c^^^A^u 
are replaced analogously. 

Note that (8) can be factorized exactly as indicated in (4): 


u?"^^ = (1-0 A )(l-0 A )(l-0 A ) u? 
ijk zz yy XX ijk 


( 10 ) 


While this property in general is lost when extending the scheme to a linear or 
nonlinear hyperbolic system, it nevertheless shows that dimensional splitting 
basically is a sound idea. 


3. SECOND-ORDER SCHEMES 

Second-order schemes for Eq. (3) may be derived, again, by extending schemes for 
Eq- (5). In order to achieve second-order accuracy in time and space we must use a 
piecewise linear initial-value distribution (see [7]), i.e.. 


n / ^ n . 

u (x,y,z) = u^j^ + 


x-x . y-y . 

1 ^ n . 1 r n 

0 u. H ;; 0 u. + 

Ax X ijk Ay y ijk 


z-z 

k . n 
Az z^ijk* 
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’ \-h ^ ^ ’ 


(11) 


here and 6^u represent locally averaged differences of u, e.g. 

6 = ^(A + V ). 

X XX 


Combining (11) with the exact solution (7) of (5) reveals that the second-order 
scheme contains terms up to o[(At)^] when restricted to one dimension, up to 
0[(At) in two dimensions, and up to o[(At)^] in three dimensions. The two- and 
three-dimensional schemes therefore can not be equivalent to products of steps with 
the one-dimensional scheme. 


For simplicity, consider first the two-dimensional scheme; with the terms arranged 
as in (8) it reads 




- a 


A {l + h(l - a )6 - ^0 A [l + a )6 + (^ - 4 a )6 lllu^., , 

v^ y y X x*- 3 V V ^ q vjjj 


y y^ 


X X 


3 \y / N/ 

y y 


0 < 1 , 0 < 1,0 < 1 . 

X y z 


3 ^x^ xJJJ^ijk^ 

( 12 ) 


When changing this scalar scheme into a scheme for a nonlinear system, we must 
replace by 


At 

THT 


{ f[(l+^ 26 x)uijkl “ f[(l-^*x)'^ijk] } » 


(13) 


a central difference of fluxes inside of volume (i jk) . The resulting method is 
unattractive because of the many intermediate steps needed for a full update; 
moreover, it is not clear that all steps are relevant for a non-diagonalizable 
hyperbolic system. 


A slight simplification results when the third-order terms on the right-hand side 
of (12) are taken together, yielding 


0 0 A A 
X y x y 



(14) 


and then redistributed in a more convenient fashion over the expressions in braces, 
e.g. with weight h in each expression. The most drastic simplification results by 
omitting all third-order terms. This does not appear to affect the stability 
condition of the scheme. Its short-wave stability is still dictated by the first- 
order scheme (8), while its long-wave stability is slightly improved (dropping (14) 
causes the long waves to slow down). 


The three-dimensional scalar linear version of this simplemost scheme is 
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= fl-aA{l + ^(l-o)6 -WA(l-iaA)-^aA(l-^aA)} 
ijk ^ X X X y y 3 z z z z 3 y y J 


- a A {l + ^(1 -0)6 -^aA(l-^aA)-%aA(l-^0A)} 

y y y z z 3 X X x x 3 z z ^ 


z z 


- 0 A |l + ^(1 - 0 )6 - w A (1 - 4 A ) - ^0 A (1 - i 0 A )}]u?., , 

z z XX 3 y y y y 3 x x ■* -* ijk* 

0 <1, 0 <1,0 <1. (15) 

X y z 


It turns out that a two-dimensional version of this scheme for the Euler equations 
was derived independently by P. Collela (private communication, 1983) and has 
already been applied successfully to an aerodynamics flow problem in a curvilinear 
grid [8]. 


4. CONCLUSIONS 

In the preceding sections it has been shown that fully three-dimensional difference 
schemes for hyperbolic systems of conservation laws can be constructed on the basis 
of a simple convection principle* 

The resulting schemes are explicit, upwind biased and stable under a combination of 
one-dimensional Courant-Friedrichs-Lewy conditions. The latter property they share 
with dimensionally factorized schemes. 

A two-dimensional second-order scheme of this kind has been put into practice by 
Eidelman, Colella and Shreeve [S]. A numerical comparison between factorized and 
non-f actorized schemes remains to be made. 
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